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I. Introduction 

One of 0,e technology goals of NASA for 

highly efficient propulsion systems to reduce the xost j y (RLV) in 

Developments of rockets for the secon § en ^ versions of conventional engines, 

the past several years have been focuse on io- Tr ortation Program to build a crew 
However, recent changes in e n eg Soace Shuttle fleet might suggest that air- 

transportation vehicle to extend the nmmilsion system for the third generation 

breathing rockets could reemerge as a tank exceeds 

RLV to replace the Space Shuttle after . ^ , launch while the payload is only 

Sme"nte I fffer clear advantages by making vehicles lighter and more 
efficient. 

Expenence in the National Aerospace Plane Pto^ the lmej98^ 

scramjet engines can achieve high speci ic Cycle (RBCC) or Turbine Based 

Whether taking a form of Rocket Based Comb ^y } K f air-breathing 
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mixing layer decreases as the Mach number increases. 

A factored form of the Gauss-Seidel relaxation 

hypersonic flow research smce its first apphcabon t ^non- qudi bnrnn^ 
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IL Numerical Method 

Let t be time; O the vector of conserved variables; E, F, and G the convective flux 
vectors; E v , F„ and G v the flux vectors for the viscous terms. The source term S 
represents production or destruction of species due to chemical reactions^ Then the three- 
dimensional Navier-Stokes and species transport equations in generalized curvilinear 
coordinates (£,, p, q) can be written as 


d,Q + dt(E - E v ) + 8 t{F- F v ) + ddG -G v )-S 


( 1 ) 


The governing equations are integrated in time for both steady and unsteady flow 

calculations For a steady-steady solution, the use of a large time step leads to a a 

convergence. For a time-accurate solution, it is desirable that the time step is detemmed 
by the physics rather than the numerical stability. An unfactored implicit scheme c 
obtained from a nonlinear implicit scheme by linearizing the flux vectors about e 
previous time step and dropping terms of second and higher orders. 


[I + aAt(DiA + DnB + DX- H)]SQ = RHS 


( 2 ) 


where 


RHS = -At[Ds(E - Ev) + D,(F - Fv) + Df{G - Gv) - S ] 

I is the identity matrix and SO denotes the correction. A, B, C, and H are the Jacobian 
matrices of the convective flux vectors and the source tenn respectively. Artificia 
dissipation models augment finite volume formulation of the right hand side. Direct 
inversion of a large block banded matrix becomes impractical m three dimensions 
because of the rapid increase of computational work and the large storage requirement. 
The LU-SGS scheme is one of the approximate factorization methods to alleviate the 

difficulties in three dimensions. 


( 3 ) 


LD~ x JJSO = RHS 


( 4 ) 


where 


L = 1 + aAt(D;A + + D r B + + DX + - A~ -B~ -C~ -H) 
D = J + aAt(A + + B* +C* - A~ -B~ -C~) 

JJ = I + aAt(D‘ ; A~ +D~B~ +D*C~ + A* + B* +C + ) 


( 5 ) 


Recent discussions 5 suggest that a different treatment of the source Jacobians might 
increase numerical stability. For example, 


L = I + aAt(D~A + +D;B* +D;<T - A~ - B~ -C H) 
D=I +aAt(A* + B* +C + -A" - ET -C -“#) 

U = / + aht{D\A +D' n B' + D*C~ + A" +B + +C -”#) 


Here, the split Jacobian matrices are constructed approximately using spectral radii to 
yield diagonal dominance. 

In the thermodynamic model, mass-averaged Cp is assumed to be a fourth order 
polynonhal 1 function of temperature with the coefficients obtained from curve fit data. 

GaJ mixture’s Cp is obtained from mass concentration weighting^ Internal energy oft 
gi i^xture is determined using mixture properties 

using the reference enthalpy (the heat of formation of species) and the integral of Cp with 
respect to temperature. Then the temperature is found from the internal energy using 
Newton iteration since Cp is assumed to be a function of temperature only. Pressure is 
obtained from Dalton’s law of partial pressure. Equivalent specific heat ratio (sensible 
enthalpy/intemal energy) is used for the Jacobian matrices. 

In the transport property model, viscosity and thermal conductivity for individual species 
,aS usingihe fourth order polynonual expressions Viscosity^ and eonducr.v.ty 
for the mixture are obtained from Wilke’ s mixing rule. Binary mass diffusivny is 
modeled usmg Chapman-Enskog formula for dilute gases m conjunction with l Lennard- 
Jones intermolecular potential energy functions. The inter-collision parameter th 
effective collision integral factor, and the effective temperature are approximated. The 
diffusion velocity of each species is calculated usmg Fick’s law and the approximated 
coefficients for binary diffusion between each species and the mixture. 

In the finite -rate-chemistry model, the rate of change of species concentration is defined 
by the movement of each species summed over all of the contributing reactions (the 
of mass action) Forward and backward reaction rate coefficients are calculated from the 
Airt^us fom. A nine-species eighteen-step reaction nu>del for 

Influence of other bodies on a particular reaction is modeled usmg a thud body efficie y 
correction. 


III. Preliminary Results 

In the first test case, hydrogen fuel is injected transversely to incoming supersonic flow 
of air The incoming air speed and temperature are Mach 2 and 1,000 K. Hydrogen is 
injected at the sonic speed while non-catalytic wall is cooled at 600 K. In Figure a and 
lb temperature contours for non-reacting and reacting flows are compared. Increased 
pressure near the exit due to combustion in Figure 2b can be compared to pressure 
contours for non-reacting flow in Figure 2a. 

When the air speed is increased to Mach 4 with 1,300° K, a lot of fuel exits without 
combustion due to insufficient mixing as shown in Figure 3a. It is shown m Fibres 3b- 
3d that various multiple injectors might enhance combustion efficiency. Hie ettec ot 
different implicit schemes on the robustness and efficiency of the reacting flow so ver or 
supersonic combustion problems will be discussed in the final manuscript. 
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Figure 2a. Pressure for frozen flow at 
Mach 2 


Figure 2b. Pressure for reacting flow at 
Mach 2 



Figure 3d. Temperature for parallel-serial 
Figure 3c. Temperature for parallel injectors at Mach 4 

injectors at Mach 4 
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